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ABSTRACT 

State  space  models  describing  the  energy  spectrum  of 
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TNTRnmirnnN 

Battery  alignment  is  one  o-f  the  critical  -factors  in 
the  control  o-f  the  weapons  systems  on  board  naval  warships. 
The  spatial  relations  o-f  the  guns,  directors,  and  launchers 
are  determined  by  reference  to  lines  and  planes  -fixed  with 
relation  to  the  hull  and  the  determination  o-f  these 
relations  constitutes  the  battery  alignment.  This  alignment 
only  compensates  -for  the  static  relations  between  the 
system  elements,  and  so,  when  these  planes  and  lines 
distort  in  a  seaway,  the  alignment  is  no  longer  correct. 

Normally,  the  reference  lines  used  -for  the  alignment 
are  the  ship  centerline  -for  train  and  the  plane  o-f  the 
primary  element,  usually  the  director,  -for  elevation.  When 
the  alignment  is  done,  the  ship  is  alongside  a  pier  in  a 
known  condition  o-f  loading  and  in  calm  seas.  At  that  time 
the  relations,  vertical,  horizontal,  and  transverse, 
between  the  directors  and  the  weapons  are  measured  and 
recorded,  and  the  -fire  control  system  is  adjusted  so  it 
compensates  -for  static  parallax  when  making  the  -fire 
control  solution.  Once  the  ship  goes  to  sea,  these 
relations  change  as  -fuel  is  burned,  stores  used  and  the 
ship  moves  in  waves.  The  sea,  in  addition  to  causing  rigid 
body  motions  o-f  the  ship  hull  in  heave,  pitch,  surge,  sway, 
roll  and  yaw,  also  causes  it  to  -flex  in  the  vertical  and 
horizontal  planes  and  to  twist  about  its  longitudinal  axis. 
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This  -flexing  distorts  the  re-ferance  lines  established  in 
the  battery  alignment  and  introduces  errors  in  the 
targeting  solution  arrived  at  by  the  -fire  control  system. 
These  errors  are  due  to  dynamic  parallax,  resulting  -from 
relative  target  location  due  to  ship  attitude  and  -flexure. 

Studies  by  the  Naval  Sea  Systems  Command,  under 
contract  to  Rockwell  International  C13,  have  shown  that  the 
■flexure  o-f  the  ship  hull  o-f  monohull  vessels  appears  to  be 
si  gn  i  -f  i  can  t ,  relative  to  rigid  body  motion.  It  seems  that  a 
method  to  measure  and  compensate  -for  these  errors  is 
necessary  to  improve  the  -fire  control  solution  and  so  to 
increase  the  -fighting  e-f-f  ec  t  i  veness  o-f  the  ship. 

Measurement  o-f  the  -flexure  could  be  done  directly,  -for 
example  by  using  a  system  o-f  sighting  instruments,  and  the 
results  -fed  to  the  -fire  control  system  as  they  happen. 
Suitable  instruments  might  be  Lasers,  light  beams  or  strain 
gauges.  All  o-f  these  require  the  primary  equipment,  along 
with  its  support  systems,  which  are  subject  to  -failure  and 
in  addition  the  system  must  be  aligned  and  kept  in 
alignment.  This  results  in  an  increased  maintenance  load  on 
ship's  personel,  and  the  reliability  problems  inherent  in 
mechanical  devices.  In  addition,  this  method  does  not 
provide  an  easy  capability  to  predict  what  the  -flexure 
might  be  at  some  -future  time,  thus  allowing  correction  o-f 
the  -fire  control  solution  at  the  actual  instant  o-f  -fire. 

An   alternative   method   involves   use   o-f  a  -filter  to 
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extract  de-flection  i  nf  ormat  i  on  -from  parameters  already 
measured  on  board  the  ship.  The  ship's  gyroscope  provides 
ship  heading  and  speed  in-formation  and  in  addition  it  can 
provide  heave,  pitch,  roll,  sway,  and  yaw  in-formation. 
In-formation  regarding  wind  speed  and  direction  is  also 
available.  Additional  in-formation,  i -f  needed,  might  be 
-found  -from  accel  erome  ters  placed  at  the  bow,  the  stern  or 
other  critical  locations.  The  advantage  o-f  using  solid 
state  devices,  instead  o-f  directly  measuring  the 
de-flections,  is  increased  reliability  and  descreased 
maintainace.  In  addition,  the  -filter  -formulation  lends 
itsel-f  to  prediction  o-f  -future  distortions,  which  can  be 
■fed  to  the  -fire  control  system. 

The  objective  o-f  this  project  is  the  -formulation  o-f  a 
mathematical  model  o-f  the  ship/sea  system  that  can  be  used 
with  a  -filter,  the  -formulation  o-f  a  -filter  to  estimate  the 
de-flections,  and  a  predictor  to  estimate  -future 
distortions.  The  equations  discribing  the  ship  can  be 
written  as  a  set  o-f  -first  order  di -f -f  eren  t  i  al  equations, 
which  lend  themselves  to  a  state-space  model  formulation. 
It  is  then  natural  to  use  a  Kalman  -filter  for  the 
estimator.  This  method  has  already  been  used  to  estimate 
the  rigid-body  motions  o-f  a  ship  for  input  to  an  automatic 
landing  system  [2,3,43. 
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CHAPTER  1;  FORMULATION 

1.1  Literature  Overview. 

Uave  induced  ship  vibration  is  the  subject  o-f  current 
interest  to  many  investigators.  The  beginnings  o-f  the 
majority  o-f  the  recent  work  dates  -from  the  late  1968's, 
when  the  per-formance  o-f  large  bulk  carriers  was 
investigated  by  Aertssen  [53,  Aertssen  and  De  Lembre  161, 
Bell  and  Taylor  C73,  and  Goodman  CS].  These  papers  were 
mainly  concerned  with  the  statistical  nature  o-f  the 
stresses  and  bending  moments  and  their  prediction, 
Aertssen,  and  Bell  and  Taylor  doing  -full  scale  mesurements 
and  Goodman  taking  a  theoretical  approach. 

The  general  trend  in  the  theoretical  approaches  has 
been  to  use  a  modal  decomposition,  as  done  -for  example  by 
Bishop  and  Price,  et.  al .  [9,18,11,12,133,  Chen  [143, 
Jensen  and  Pedersen  [153,  and  Kagawa  and  Fuj i ta  [163.  This 
has  an  advantage  when  working  with  regular  waves  and  for 
seperating  out  the  contributions  of  the  various  modes. 

The  problem  o-f  the  exciting  hydrodynamic  -forces  has 
been  addressed  by  Salvesen,  Tuck,  and  Faltinsen  [173  -for 
the  rigid  body  motions,  using  strip  theory,  providing 
linear  expressions  -for  the  distribution  o-f  the  wave  -forces 
over  the  hull.  Skj8rdal  and  Faltinsen  [183  have  developed  a 
linear  theory  -for  the  hydrodynamics  o-f  a  springing  hull. 
Jensen   and  Pedersen  [153  provide  a  non-linear  -formulation. 
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Betts,  Bishop  and  Price  CI?]  tie  the  strip  theory   approach 
together  with  structural  representations. 

The  use  o-f  state  space  techniques  -for  the  prediction 
o-f  hull  vibration  does  not  seem  to  have  been  used  by  any 
researchers  as  no  references  were  -found  by  the  author 
discribing  this  approach. 

1.2  Ship  Model 

There  is  no  establ ished  procedure  on  how  to  best  model 
the  ship  and  the  sea  in  order  to  predict  flexural 
distortions.  One  major  consideration  is  the  size  of  the 
models  which  can  make  the  use  o-f  state-space  techniques 
very  cumbersome.  For  this  reason,  -for  example,  it  is  best 
not  to  use  exact  structure  and  hydrodynamic  models,  such  as 
■finite  element  techniques,  and  instead  to  try  to  reach  a 
compromise  between  model  size  and  modeling  accuracy.  A 
question  arises  when  trying  to  combine  a  structural  model, 
such  as  modal  representation,  with  a  random  sea.  The  sea, 
because  o-f  dispersion,  has  a  wavelength  which  is  a  function 
o-f  wave  -frequency.  Even  though  the  random  sea  may  be  easily 
represented  by  Fourier  techniques  at  a  speci-fic  point,  its 
representation  at  some  other  point  in  space  would  be  the 
result  o-f  the  propogation  o-f  an  in-finite  number  o-f 
■frequency  components.  How  to  include  these  phase  shi-ft 
e-f-fects  in  a  modal  analysis  is  not  clear.  In  addition,  the 
computational   load  may  be  large  due  to  the  number  of  modes 
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required  to  accurately  represent  the  quasi-static  hull 
distortion,  given  that  the  natural  -frequencies  o-f  the  hull 
are  much  higher  than  the  wave  -frequencies  containing  the 
major  part  o-f  the  wave  energy. 

The  approach  used  here  is  to  approximate  the  hull  as  a 
quasistatic  non-uniform  beam.  Using  simple  beam  theory,  one 
starts  by  considering  a  short  section,  dx ,  o-f  the  hull  (see 
■figure  1.1).  The  governing  equation  o-f  the  beam  vibrating 
in  a  plane  o-f  symmetry  is  -found  by  summing  the  vertical 
■forces  acting  at  point  "a" 

m(x)  dx  d/dt<dy/dt)  =  -f<x,t)  dx  -  V<x,t)  +  U(x,t)+dV/dx  dx 

=  -f(x,t)  dx  +  dV/dx  dx  (1.1) 


or 


m  d/dt<dy/dt)  ■  -f  +  V  <1.2) 


Similary,  summing  the  moments  about  "a"  one  obtains 


M'<x,t)  =  -V<x,t)  <1 .3) 


Then    by    using    (1.22)     in    (1.21)    one    -finds 


m    d/dt(dy/dt)    =    -f    -  M"  <1.4) 
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F(Xjt)  a  x 


Fi  gure  1.1 
Beam  Element 


-11- 


From  simple,  linear  beam  theory,  one  has  the  relation 


M  -  Ely"  (1.5) 


Using  (1.24)  in  (1.23)  one  obtains 


m  d/dt(dy/dt)  =  i    -  (Ely")"  (1.6) 


which   is  the  governing  equation  -for  vertical  vibrations  o-f 
a  beam . 

I -f  damping  is  allowed,  then  one  has 

m  d/dt(dy/dt)  +  b  dy/dt  +  (Ely")"  =  -f  (1.7) 


Once  one  has  the  governing  equation,  the  next  step  is 
to  express  it  in  a  state-space  -form.  To  do  this,  let  xl  =  y 
and  x2  =  dy/dt.  Then  with  this  coordinate  trans-formation, 
( 1 .26)  becomes 


dxl/dt   =   x2 

(1  .8) 
dx2/dt   =    -(EIxl")"    /m-bx2/m+-f 


The  next  step  is  to  represent   the   space   derivative. 
Using   a   -finite   di-f-ference  approach,  one  lets  h  =  dx ,  the 
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distance  between  points  on  the  hull.  Then  i  -f   there   are   N 
points, 

L  =  (N-l)h  (1 .9) 


where   L   is   the   length   o-f  the  ship.  Expanding  the  space 
derivative,  one  has 

(EI  xl">"  =  EI"  xl"  +  2  EI'  xl'"  +  EI  xl""    (1.18) 


Now  each  part   can   be   represented   seperately   by   -finite 
di -f-f  erences,  with  EI,  EI',  and  EI''  as  input  data. 

At   the  ends  o-f  the  ship,  the  conditions  o-f  zero  shear 
and  zero  bending  moment  lead  to 


EI  xl"  =  8 

(1.11) 

(EI  xl")'  =  8 


=  EI'  xl"  +  EI  xl'" 


Equations  (1.29)  with  boundary  conditions  (1.291)  can 
be  represented  to  the  order  desired  and  will  result  in  the 
■formation  o-f  a  set  o-f  equations,  written  in  matrix  -form  as 
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dXl/dt    =   A1*X1     +    B1*E 
Y        =    C1*X1 


<1 .12) 


where 


X'    =    <xll ,xl2,x21 ,x22, . . . ,xNl ,xN2) 

xij    =    x    sub   j    at    point    i 

Y   =    vector    o-f    responses    at    the    N    points. 

A   =    a    2N    x    2N   matrix    containing    the    -finite 
di-f-ference    representation. 

E'    =    (-F1  ,-f2,  .  .  .  ,-fN) 

B  =  distributes  the  -fj's  to  the  proper  state 
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1 .3  Sea  Spec trum 

Ocean  waves  are  a  function  o-f  both  the  wind  speed  and 
duration,  among  other  parameters.  One  o-f  the  commonly  used 
wave  spectra  that  attempt  to  account  -for  these  two  -factors 
is  the  Bret Schneider  Spectrum,  in  which  the  intensity  is 
speci-fied  by  the  significant  wave  height,  H,  and  the 
duration  by  the  modal  wave  -frequency,  Wm.  This  spectrum  has 
the  -followinQ  -form. 


S<w)  =  1.25  H**2  wm**4  exp  i-1  .  25<U)m/w)  **4>      (1.13) 
4         w**5 


The  movement  o-f  a  ship  through  the  water  causes  it  to 
experience  excitation  not  at  te  wave  -frequency,  w,  but  at 
the  encounter  -frequency,  we,  where 


we  =  w  -  kU  cos<beta>  <1.14) 


The   energy  o-f  the  waves  remains  the  same,  however,  so 
that 
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S<w)dw  =     S<we)dwe 

(1 .15) 
—    S<w)  dw   dwe 
dwe 


there-fore 


S(we)  =    S(w) 


dwe  <1 .16) 

dw 


where  the  modulus  o-f  dwe/dw  is  used  because  S(we)  is  single 
si  ded. 

Tr  i  an  ta-f  y  1  1  ou  and  Athans  C23   have   -found   a   rational 
approximation  to  (1.38)  in  the  -following  -form 


S(we)  =  1  .25  H**2  B<alpha) we**4 (1.17) 

4Wm  <  l+<we/wo)**4)**3 


where 


-16- 


alpha  =  U_wm  cosCbeta}  (1.18) 

9 

B<alpha)  =    1 .9339  (1.19) 

1+2  alpha 

wo(alpha)  =  wm(l+alpha)  (1.29) 

9.8489 


The  corresponding  causal  system  is 


H<s)  =   So  s**2 (1.21) 

( l+2J(s/wo)+(s/wo)**2)**3 


where 


So  =  1 .25  H**2  B(alpha)  (1.22) 

4     Wm 

J  =  8.787 


By  writing   (1.38)   in  standard  state  space  -form  (see 
chap  2)  one  has 


dXs/dt  =  As*Xs  +  Bs*p 

(1 .23) 
ys  =  Cs*Xs 
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where  p  is  white  noise  o-f  intensity  q  and 


As  =  : 


Bs'  = 

Cs  = 


I       9 

1 

8 

9 

8 

9 

:-wo**2 

-2Jwo 

9 

wo**2 

8 

8 

9 

9 

9 

1 

9 

8 

<1 .24a) 

9 

9 

-wo**2 

-2Jwo 

9 

wo**2  : 

9 

9 

9 

9 

9 

1 

8 

9 

9 

9 

-wo**2 

-2Jwo  : 

[   9 

9 

9 

9 

9   p 

i  *wo**3: 

1    <1.24b) 

C   So 

9 

9 

9 

9 

8 

1    <1.24c) 

Since  white  noise  has   a   constant   spectrum   tor   all 
f requenc i  es 


Sp(w)  =  q 
P" 


<1 .25) 


and 


Sp(w)  =  H<w)  **2  S<w) 


=  wo**3  So 


we**4  wo**4 


<1 .26) 


(wo**4  +  we**4)**3 


one  may  ideniti-fy  the  intensity  as 


q  =  p  i  wo**3  So 


(1 .27) 
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and  the  driving  noise  is  taken  as 


p  -  N<8,1) 


<1 .23) 


with  pi  wo**3  So  included  in  the  model  B  matrix. 
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1 .4  Sea  Model 

The  representation  o-f  the  sea  must  take  account  o-f  the 
dispersion  ot  the  waves  in  space  and  still  allow  -for  their 
random  nature.  I*  the  relation  between  driving  point  1,  pi, 
and  distant  point  2,  p2,  is  dependent  on  a  -fixed  wave 
number,  then  the  sea  at  p2  will  not  be  properly  related  to 
pi.  To  -find  a  relation  that  will  give  the  correct  relation 
between  pi  and  p2,  one  starts  with  Laplace's  equation  and 
the  -free  surface  condition: 


Qxx  +  Qzz  =  9  in  the  -fluid 

<1 .2?) 
Qtt  +  g  Qz  =  8  on  the  -free  sur-face 


Then  by  assuming  that  the  potential  can  be  written  as 


Q(x,z,t)  =  U<x,t)  exp<kz>  <1.38> 


and  substituting  (1.41)  into  the   -free   sur-face   condition, 
one  -finds 

Qtt  +  gQz  =  <Utt  +  gkU)  exp<kz>  (1.31) 


I-f  the  deep  water  dispersion  relation,  k  =  w**2/g,  is  used, 
(1.42)  reduces  to 
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Utt  +  w**2  U  =  8 
or 

Utt  -  Utt  =  0 


showing   that   the   -free   sur-face  condition  is  satis-fied  as 

expected.   Then   substituting   < 1 . 41 >   into  the    Laplace 
equation,  one  obtains 

g**2  Uxx  +  Utttt  =  8  (1.32) 


It  is  intresting  to  note  that  (1.43)  may  be  obtained 
■from  the  wave  equation  by  using  the  expression  -for  wave 
phase  velocity  in  deep  water. 

c  =  w  /  k 
and 

Utt  -  c**2  Uxx  =  8 
-from  wh  i  ch 

Utt  -  w**2/(w**2/g)**2  Uxx  =  8 
or 

w**2  Utt  -  g**2  Uxx  =  8 
and  since  w**2  U  =  -d/dt(dU/dt) 

Utttt  +  g**2  Uxx  =  8 

Equation       (1.43)       can      now      be      represented      by   -finite 
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di-f-ference  equations  in  space  and  as  a  set  o-f   -first   order 
equations  in  time  as 


dXv/dt  ■  Av*Xv  +  Bv*ys 
E    =  Cv*Xv  +  Dv*ys 


(1 .33) 


where 


Xv'    -    <ull ful2,ul3,ul4,u21 , . . . ,uN4) 

u j 1    =   Uj    at    the    j-th    point 

uj2  =    dUj/dt 

uj3  =    d/dt<dUj/dt) 

uj4  =   d/dt<d/dt<dUj/dt>> 

E   =    vector    o-f    wave    elevations    at    the   N   points 
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1 .5  Force  Model 

The  hydrodynamic  exciting  -forces  on  a  ship  hull  by  sea 
waves  have  been  studied  -for  many  years.  The  approach  that 
seems  the  most  -fruit-ful  -for  this  project  is  the  strip 
theory.  0-f  the  several  -forms  available  116,191,  the  method 
o-f  Gerritsma  and  Beukelman  C193  will  be  used  here. 

The  hydrodynamic  exciting  -force  on  the  ship  may  be 
di  scribed  as  being  the  result  o-f  the  di-f-ference  between  the 
de-flection  o-f  the  ship  hull  and  the  wave  pro-file,  z=y-w. 
Then  Gerritsma  and  Beukelman  have  proposed  that  the  -force, 
F<x,t),  acting  on  the  ship  may  be  written  as 


F(x,  t)  =  - 


D 
Dt 


Ma(x) 


Dz 
Dt 


N<x)  *  _Dz_  +  rgB<x)z 
Dt 


where 


D  =  d   -  d   dx  =  d  r  -  U  d 
dt    dx  dt    dt      dx 


Then,  using  z=y-w ,  and  writing  the  -force  as 


<1 .34) 


where 

H<x,t) 


F<x,  t)  =  -H<x,t)  +  Z<x, t) 


Ma(x)  D_2   +   N<x)-U  dMa   Dy  +  rgB(x)y 
DtA  dx    Dt 


Z<x,t)  =  Ma(x)  Daw   +   N<x)-U  dMa   Dw  +  rgB(x)w 

Dt*-  dx    Dt 


U  .35) 


<1 .36) 


a  .3?) 


the  force  can  be  seperated  into  the  part  due   to   the   hull 
motions  and  the  part  due  to  the  waves.  The  part  due  to  the 
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hull  vibrations  may  be  moved  to  the  le-ft  hand  side  o-f  the 
equation  and  included  in  the  structural  model.  The  wave 
exciting  -forces  then  may  be  written  in  a  -form  that  may  be 
used  in  the  state-space  -formulation.  Then  the  C-matrix  and 
D-matrix  o-f  the  sea  model  are  modi-fied  with  equation  (1.53) 
to  output  the  -force  instead  o-f  the  elevation.  The  output 
matrix  o-f  the  spectral  model  is  modi-fied  to  output  the 
•first  and  second  time  derivatives  o-f  the  wave  elevation  in 
addition  to  the  elevation  itsel-f. 

1 .6   Total  model 


The  models  -for  the  sea  spectrum,  the  sea  surface,  and 
the  ship  can  be  combined  into  one  model  by  identifying  the 
output  o-f  one  with  the  input  o-f  another  in  cascade  -fashion. 


where 


dX/dt  =  A*X  +  B*p 
Y  =  C*X 


A      - 

:  ai 
:   e 
:  0 

Bl*Cv 
Av 
8 

Bl*Dv*Cs    : 
Bv*Cs 
As 

B'    = 

[      0 

8 

Bs'          ] 

C      = 

[    Cs 

8 

0                 ] 

X'    = 

(    XI  ' 

Xv' 

Xs'            ) 

<1 .38) 


p  ~  N(0,1) 
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CHAPTER  2:   STATE  SPACE  TECHNIQUES 

2.1  State  Space  Models 

In  order  to  apply  Kalman  Filtering  techniques,  one 
must  -formulate  the  problem  in  the  -form  o-f  a  state-space 
model.  This  is  a  set  o-f  -first  order  di -f-f eren t  i  al  equations 
that  discribe  the  dynamics  o-f  the  system,  its  inputs,  and 
i  ts  ou tpu  ts. 

The  state  o-f  a  system  may  be  de-fined  at  any  time  as  a 
minimum  (not  neccessarily  unique)  set  o-f  parameters  which 
must  be  specified  to  discribe  that  system.  A  state  vector 
is  a  vector  composed  o-f  these  parameters  (state  variables). 
The  number  o-f  state  variables  required  is  equal  to  the 
order  o-f  the  system,  n. 

For  example,  a  simple  mass-spring-damper  system  is 
discribed  by  a  second  order  d  i -f-f  eren  t  i  al  equation. 

m*d/dt(dz( t)/dt)  +  b*dz(t)/dt  +  k*z(t)  =  -f(t)    (2.1) 

Then  one  may  de-fine  the  state  variables  as  the 
position,  z,  and  the  velocity,  dz/dt,  o-f  the  mass.  Next  let 
xl  =  z  and  x2  =  dz/dt.  With  this  change  o-f  variables,  (2.1) 
may  be  wr i  t ten  as 


dxl/dt   =   x2 

(2.2) 

dx2/dt    =    -k/m*xl    -    b/m*x2    +    -f/m 
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or  in  matrix  notation 


dX/dt  =  A*X  +  B*F 


(2.3) 


where 


A  ■ 

B'    - 
X'  = 


8     1 
-k/m  -b/m 
C    8     1   ] 
<   xl    x2   ) 


Alternatively,  because  the   state   variables   are   not 

unique,   one  may  define  the  state  variables  as  the  -force  in 

the   spring,   K,   and   the  momuntum   o-f   the   mass,   P.   A 

trans-formation  may  be  -found  -from  one  to  state  vector  to  the 
other  by  using 


K  =  k*z 

P  =  m*dz/dt 


<2.4) 


One  may  also  de-fine  an  input  vector  U,  containing  the 
r  inputs  to  the  system  and  an  output  vector,  Y,  containing 
the  m  outputs.  The  system  equations  are  then  a  set  o-f  n 
•first  order  di-f-ferential  equations 


dxl/dt  =  -fl(X,U,t) 


dxn/dt  =  -fn<X,U,t) 


(2.5) 
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or 


dX/dt  =  F<X,U,t)  (2.6) 

The  output  equations  relate  the  output  vector,   Y,   to 
the  state  vector  and  the  input  vector 

y\    =  gl<X,U,t) 


■  ■ 


yrri  =  gm  <  X  ,  U  ,  t ) 
or 


Y  =  G(X,U,t)  <2.7) 
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I-f  the  system  is  linear  or  linearized,  one  may  write 
the  -functions  -f  i  and  gi  as  linear  combinations  o-f  the 
states  and  the  inputs. 


j*xj  +  ^/__j   bij*uj 


dxi/dt  =       ai 

J=l  j=l 


n 


yk  =  ^/        ckj*xj  +  /        dkj*uj 
3=1  J=: 

i  =  CI ,n3 
k  =  CI ,m] 


(2.8) 


(2.9) 


In  matrix  notation  this  becomes 


dX/dt  =  A*X  +  B*U 
Y    =  C*X  +  D*U 


<2.18> 


where  A  i s  an  n  x  n  matrix  discribing  the  system 
dynamics,  B  is  n  x  r  and  di  scribes  the  distribution  o-f 
inputs,  C  i  s  an  m  x  n  matrix  discribing  the  distribution  o-f 
states  to  the  outputs,  and  D  is  a  m  x  r  matrix  discribing 
the  distribution  o-f  inputs  to  outputs.  The  -form  discribed 
by  <2.18)  is  known  as  Standard  State-space  Form. 

2.2  Kalman  Fi 1 ter 


In  this  problem  the  object  is  to  be  able  to  estimate 
the  -flexural  response  o-f  the  ship  hull  at  time  t  and  -from 
this   to  predict  what  the  response  will  be  at  time  t+dt.  To 
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do  this,  one  common  method  is  to  use  the  Kalman  Filter. 
This  is  a  method  which  derives  an  optimal  <in  the  least 
squares  sense)  estimator  -for  the  state  o-f  a  system 
represented  by  a  state-space  model . 

Assume   one   has   a  mathematical   model  o-f  the  system 
under  consideration, 


dx/dt  ■  A*x  +  B*w 

<2.11) 
P   =  C*x 


where  w  is  white  noise  o-f  intensity  Q,  the  expected  value 
o-f  x  at  t=8  ,  ECx<8)],  is  x8  =  0,  and  the  covariance  o-f  x  at 
t=8  ,  ECx<8)x'  (8)  ]  ,  is  X8 .  Assume  -further  that  one  has  1 
measurements,  z,  which  are  related  to  x  by 

z  =  M*x  +  v  (2. 12) 

where  v  is  measurement  noise,  modeled  as  white  noise  of 
intensity  R,  uncorrelated  with  the  process  noise,  w. 

Also  assume  that  the  system  is  observable,  meaning 
that  the  initial  condition  o-f  any  state  may  be  -found  by 
measuring  the  output  -for  a  -finite  period  o-f  time. 

Then  one  may  propose  an  observer  o-f  the  -form 

dy/dt  =  A*y  +  K*<  z  -  M*y  )  (2.13) 

y(8)  =  EC  x<8)  ]  =  8 

K  is  known  as  the  Kalman  gain  matrix. 

-29- 


One  may  then  de-fine  the  error  in  the  state  estimate  as 


e  =  y  -    x 
and  the  error  covariance  as 
P  =  EC  ee'  ] 


(2.14> 


(2.15) 


It   may   be   shown   [23]   that   the  equation  which  controls 
propagation  o-f  the  error  covariance  is 


dP/dt  =  <A-KM)P  +  P(A-KM)'  +  KRK'  +  BOB' 
P(0)  =  X(S)  =X9 


(2.16) 


Examining  the  covariance  matrix,   P,   -finds   that   the 
diagonal  is  made  up  o-f  the  variances  o-f  each  o-f  the  states, 


P  =  EC  ee'  3  = 


el 
e2 


en 


C  el  e2  ...  en  ] 


el**2  el*e2 
e2*el  e2**2 


en**2 


(2.17) 


So  i -f  one  takes  he  Trace  o-f  P 
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n 

s 

i  =  l 


■2 


one  has  the  sum  o-f  the  squares  o-f  the  errors  in  each  o-f  the 

states.    Athans   C233  shews,   using   the   matrix   minimum 

prinicipal,  that  the  K  that  minimizes  the  Trace  o-f   P,   ie; 

the  sum  o-f  the  squares  o-f  the  error,  is 

K  =  PM'R  (2.19) 

I-f,  as  in  the  present  case,  none  o-f  the  matrices  A,B,M,Q 
and  R  are  -functions  o-f  time,  then  dP/dt  — >  9  i  n  the  steady 
state  solution,  and  one  has 

9  =  AP  +  PA'  +  BQB'  -  PM'R  MP  <2.28> 

Summarizing,    then,    i-f   one   has   a   linear   system 
descr  i  bed  by 

dx/dt  =  A*x  +  B*w 

<2.21> 

p    =  Cx 

one  may  estimate,  with  minimum  square  error,  the  state  x  by 
using  y  which  is  described  by  the  equation 

dy/dt  =  A*y  +  K(z-My)  <2.22) 

K   =  PM'R  <2.23> 
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with  P  given  by  (in  the  steady  state  case) 

9  =  AP  +  PA'  +  BQB'  -  PM'R  MP  (2.24) 

So-ftware  packages  exist,  for  example  CTRLC  1261  ,  which  will 
solve  2.24  -for  the  steady  state  values  in  P  and  so  -for  the 
val ues  in  K. 

Once  the  values  o-f  the  state  at  time  t8  are  known, 
either  -from  direct  measurement  or  -from  estimates,  the  next 
task  is  to  predict  the  the  states  at  some  -future  time, 
t-f=t8  +  dt.  To  do  this,  one  may  solve  2.291  subject  to 
x<t8)=x8  as 


x<t)  =  exp<A< t-t8>}*x8  + 


Nw,  suppose  that  at  time  t8  ,  one  wants  to  predict 
x< t8+dt) .  Then 

r 

x<t-f)  =  exp{A<dt)>*x8  +   \   expCA<  t-f-s)  >Bw<  s)ds     <2.26) 

(J  t8 

but  t8<s<=t-f  is  some  -future  time  and  w<s)  is  unknown 
because  w  is  random.  In  -fact,  since  w  is  white  noise  and  so 
completely   random,   it   is   completely  unknown  and  the  est 
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estimate  o-f  it  is  its  average  value  or  expected  value, 


EC  w< t)  ]  =  8 


<2.27) 


There-fore,    the    best    estimate    o-f    x(t-f)     is 


xp(t-f)    =    exp<Adt>*y<t8> 


(2.28) 


where  y    has  been  used  in  the  initial  condition   because   it 
is  obtained  -from  the  Kalman  -filter. 
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2.3  Approximations 

The  question  o-f  the  order  o-f  the  -finite  di-f-ference 
representation  will  be  addressed  next.  The  accurarcy  o-f  the 
distortion  results  are  dependent  on  that  order.  For  just 
one  point  o-f  the  ship  model,  say  j,  equation  1.27  can  be 
written  as  (with  j  understood  in  the  coe-f-f  i  c  i  ents) 


dxjl/dt 
dxj2/dt 


=  xj2 


<2.29) 


-EI  xjl'  '" 
m 


-  EI  '  xjl"'  - 


EI"  xjl"  - 


m 


m 

b  xj2  +  jF^ 
m        m 


Since  the  highest  derivative  in  x  is  the  controlling  -factor 
in  the  discretation  error,  the  case  o-f  a  uni-form  ship,  with 
EI'  and  EI''  both  equal  to  zero  will  be   investigated.   The 
method  used  will  be  that  used  by  Smith  C233. 
At  point  j,  (2.38)  becomes 


dX/dt  =  A*X""  +  D*X  +  F 


<2.38) 


where  the  matrix  A  is 

8       8 
:  -El/m     8 


A  = 


and  D  i  s 


D  = 


8       8 
8     -b/m 
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and    F    i  s 


F   = 


8 
■f/m 


Now  write    the    -fourth    derivative    o-f    X    as 


X<j>""    ■   XCj+2)    -4X<j+l)    +    6X<j)    -    4X<j-l)    +    X(j-2) 

h**4  (2.31) 


Using    <2.32)     in    (2.31),    one    -finds 


dX(j)/dt    =    A   XCj+2)     -4AX(j-H)     +    6A   X(j) 

h**4 

-    4AXCJ-1)    +    AXCj-2) 
h**4 

+    D   X<j) 


(2.32) 


and  expanding  the  X<j+1)  by  Taylor  series 


X(j  +  1)  =  X  +  h*X'  +  J___h**2*X"  +  1  h**3*X/xx  +  .. 

2!  3! 


C2.33) 


Substituting   (2.34)   into   (2.33),   keeping  terms  to  order 
h**8  and  collecting  terms  leads  to 


dX(j)/dt  =  A*X""  +  D*X  +  F  + 

<  l*A*Xg:*h**2  )  +  (  1  *A*XZZn:*h**4  ) 
6  S3 


C2.34) 
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where  the  term  in  the  -first  bracket  can  be  interpeted  as 
the  order  h**2  error  and  that  in  the  second  bracket  as  the 
order  h**4  error. 

In  order  to  see  the  e-f-fect  on  a  particular  solution  to 

(2.38),  consider  the  solutions  o-f  the  natural   modes   o-f  a 

uni-form  beam  C223.  The  space  dependent  part  is  composed  o-f 
sinusoids  and  hyperbolic  -functions  in  the  variable  kx. 

w(x)  =  sin(kx)  +  sinh(kx)  +  s*(  cos(kx)  +  cosh(kx)  )  (2.35) 

where  s  is  a  constant  dependent  on  the  mode .  I -f  L  is  the 
length  then 

kl_  -  (j+  )*Pi 

de-fines  the  approximate  values  o-f  k.  Taking  the  eight 
dirivative  o-f  (2.35)  brings  out  a  k**3  term,  so 

w    =  k**8*w 
and  (2.31)  can  be  written  as 

dX/dt  =  A*X""  +  D*X  +  F  +  E  (2.34) 
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where 


E  -  J__  A*XZnE*h**4 
88 

=  J A*k**8*X*h**4 

88 


Then  since 


k**4*X  =  X"" 


one  may  wr i  te 


dX/dt  =  A*X""  <  1  +  _l_*h**4*k*4  )  +  D*X  +  F     <2.37) 

88 


Then  the  relation  between  the  point  spacing  and   the   error 
may  be  -found  -from  the  second  term  in  the  brackets 


e  ■  _1_  *h**4*k**4 
88 


I-f  the  maximum  desired  error  is  eo,  then 


h**4  <  88*  eo 
k**4 
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or 


h  <   2.9?*L 


eo 


(2.38) 


By   similar  means  the  relation  -for  order  h-squared  can 
be  -found  to  be 


h  <   2.45*L 


eo 


< j+f)Pi 


<2.3?) 


Table  2.1  shows  h  vs.  eo  -for  both  orders  -for  the 
springing  mode,  and  table  2.2  shows  h  vs.  j  -for  a  -fixed 
eo=8.125.  Length  is  188  meters. 

Table  2.1 
h  vs.  eo  <j  =  l) 


eo 

h**2 
h 

> 
Npts 

h**4 
h       Npts 

8.5 
8.25 
8.125 
8.8625 

36.76 
26.88 
18.36 
13.88 

3 
4 
6 
8 

53.35 
44.87 
37.73 
31  .72 

2 
3 
3 
4 
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Table  2.2 


h  vs.  j  <eo  =  9 . 125) 


1 

2 
3 
4 


h**2 


h 

Npts 

h 

Npts 

18.38 

6 

37.73 

3 

11  .83 

9 

22.64 

5 

7.88 

13 

16.17 

7 

6.13 

17 

12.58 

8 

h**4 


The  most  obvious  -fact  to  be  gained  -from  Table  2.1  is 
that  -for  a  given  accuracy,  the  number  o-f  points  and  hence 
the  number  o-f  states  is  much  less  with  the  order  h**4 
approximation.  Each  added  point  requires  six  additional 
states  in  the  total  model,  two  -for  the  ship  model  and  -four 
■for  the  sea  model.  However,  only  -five  points  are  required 
to  approximate  the  -fourth  direvative  to  order  h**2  but 
eight  points  are  required  -for  order  h**4.  These  numbers  set 
the  minimum  number  o-f  points  required  -for  any  accuracy. 

As  the  number  o-f  nodes  in  the  vibration  increases,  the 
number  o-f  points  -for  a  given  accuracy  increases  rapidly,  as 
seen  in  Table  2.2.  Again,  the  order  h**4  approximation  is 
the  less  restrictive  o-f  the  two.  It  only  requires  8  points 
to  achieve  12. 5X  accuracy  with  order  h**4,  while  17  points 
are  required  -for  the  order  h**2  approximation. 

This  means  that  while  it  is  possible  to  use  an  order 
h**2  approximation,  an  actual  i np 1 emen tat i on  will  require 
both  a  higher  order  and  a  greater  number  o-f  points. 
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CHAPTER  3;  NUMERICAL  RESULTS 

3.1  Stabi  1  i  ty  o-f  Models 

The  mathematical  models  o-f  the  sea  spectrum  and  the 
ship  structure  are  stable  as  modeled.  However  the  -force 
model  is  not.  The  -fourth  order  equation  discribing  the  sea 
surface  has  -four  solutions,  two  o-f  which  are  exponentials, 
and  it  is  these  solutions  that  can  cause  instability.  The 
equat  i  on  is 

Utttt  +  g**2  Uxx  =9  <3.1> 

which  has  solutions  in  the  -form 

U(x,t)  "    i    e  x  p  <  i  w  t )  ,  exp(-iwt),  exp<wt),  exp(-wt)  > 

<3.2> 

The  exponential  roots,  in  the  analytical  problem,  are  taken 
care  o-f  by  the  initial  conditions.  When  using  a  computer, 
however,  even  i-f  the  initial  conditions  are  set  equal  to 
zero  subsequent  round-o-f-f  errors  introduce  -fictitious 
n on -zero  initial  conditions,  which  excite  growing 
exponent  i  al s. 

To  see  this,  consider  the  simple  example 

Ytt  -  a**2  Y  =  8  (3.3) 

Y<8)  =  8 

Y<  in-f  ini  ty)  =  9 

which  has  the  analytic  solution 

Y<t)  =  A  exp<at)  +  B  exp(-at)  (3.4) 
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By    approximating    this    by   -finite    di  -f  -f  erences ,    one    obtains 


Ytt    -  Y( i+i )    -2*Y<  i )    +    Y<  i-1)    =    a**2   Y< i ) 

h**2 


or 


Y<i  +  1)  -(2  +  <ah)**2>  Y< i )  +  Y<  i -1 >  =  8 
To  solve  this  -finite  di-f-ference  equation,  take 

Y<n)  =  r**n 
then  one  obtains 

r**2  -<2  +  <ah>**2)  r  +  1  =  0 
and  then  -for  h  <  1/a 


<3.5) 


r  ~  1  +  ah  +  0<ah**2) 


and 


r  ~    1  -  ah  +  0<ah**2) 
Then  taking  a  linear  combination  o-f  the  solutions 
Y<n)  =  A  <1  +  ah)**n  +  B  <1  -  ah)**n 


(3.6) 


it  can  be  seen  that  no  matter  how  small  A  is,  i -f  it  is  not 
exactly  equal  to  zero,  the  part  o-f  the  solution  involving 
the  growing  exponenitial  will  eventually  dominate  the 
solution,  swamping  the  correct  solution.  I -f  the  objective 
were  only  to  obtain  the  solution  to  the  equation,  one  might 
try  integrating  backwards  in  time.  Here  though,  the 
objective  is  to  produce  the  solution  in  real  time,  so  the 
equation  must  be  integrated  -forward  in  time. 

This  is  not  as  severe  a  problem  as  it  might  seem  at  -first 
look.  One  o-f  the  properties  o-f  the   Kalman   Filter   is   the 
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stabilization  o-f  the  -filter  model.  The  gain  that  results 
will  result  in  a  closed  loop  observer  matrix, 

Acl  =  A  -  K*C 

that  has  all  its  eigenvalues  in  the  le-ft-hal-f  plane.  The 
only  place  the  problem  arises  is  in  simulation  o-f  the 
system.  By  de-fining  a  reconstruction  error,  e=x-xhat,  and 
using  the  expressions  -for  x  and  xhat,  it  can  be  seen  that 
the  error  must  obey  the  differential  equation 

de(t)  =  CA-KC3  e<t)  <3.7) 

dt 

This  means  that  the  reconstruction  error  goes  to  zero  as  t 
becomes  large  only  i -f  the  closed  loop  observer  matrix  is 
asymptotically  stable. 

Simulation  o-f  the  system  is  necessary  to  test  the 
observer  as  the  model  is  varied.  The  method  used  is  to  sum 
a  number  o-f  sinusoids,  each  o-f  which  has  an  amplitude 
proportional  to  the  energy  in  the  sea  spectra  at  that  wave 
■frequency  and  a  phase  that  is  random  with  respect  to  the 
phases  o-f  the  other  sinusoids. 

Assume  that  we  are  given  a  sea  spectrum,  S(w) ,  and  we 
wish  to  simulate  the  sea  elevation  at  a  number  o-f  points. 
Since  the  instantanous  sea  elevation  at  a  point  x  and 
■frequency  w  can  be  represented  as  the  square  root  o-f  the 
energy   at   that   -frequency   (multiplied  by  2  since  the  sea 
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spectra    is    only   de-fined   -from    9    to    in-finity), 

'vKxjW)    =  V  2    S<x,w>' 

then    the    elevation    at    x,     including   all       -frequencies      should 
be 


oc 


V<x)    =    ReC       \\J2    S(w)    dw '       exp{     ip<w)     >> 
8 


(3.8) 


where  p  is  the  random  phase.  Then  including  the  time  -factor 
one  has 


oo 


V<x)  =  Re<  \\J2    S<w)  dw '  expC  i  <wt  +  p<w)  >> 
8 


<3.?> 


For   a   computer   simulation,   equation   <3.9)  must  be 
discretized  to  become 


Nw 
V(x,t)  -  /       i\j  2    S<wj)  4w.j  '     cos<  wjt  +  p<wj)  >  (3.18) 
j  =  l 

where  here  p   is   a   random   number   uni-formly   distributed 
between  8  and  two  pi . 

The   simulated   sea   surface  is  then  used  to  drive  the 
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ship  model  to  -form  a  simulation  model  ,  the  output  o-f  which 
is  used  as  a  source  o-f  the  measurements  supplied  to  the 
observer  model.  This  model  was  run  and  the  results  stored 
in  a  computer  -file  -for  later  use.  Figure  3.1  shows  a 
representative  run  o-f  the  simulation. 

3.2  Estimation  and  Prediction  with  Correct  Model 

The  Kalman  -filter  -for  the  total  model  is  required  to 
construct  optimal  estimates  o-f  the  states  o-f  the  system 
that  can  be  used  to  predict  the  -future  output  o-f  the 
system.  This  includes  those  states  that  cannot  be  measured 
directly,  such  as  those  o-f  the  sea.  To  do  this  it  must  have 
a  source  o-f  in-formation  about  the  system. 

The  -filter  gains  are  a-f-fected  not  only  by  the  model 
used  but  also  by  the  noise  in  the  measurements  themselves. 
This  noise  results  -from  a  number  o-f  sources.  For  example, 
■from  measuring  instrument  noise,  slamming  o-f  the  hull, 
local  vibration  e-ffects  and  vibrations  excited  by  external 
sources.  The  -filter  then  must  reject  the  higher  -frequency 
noise  to  obtain  the  needed  in-formation  about  the  states. 
The  measurements  are  -found  -from 

z  =  y  +  v,  v  ~   N<8,R)  (3.11) 

where  R  is  the  measurement  noise  covariance  matrix,  de-fined 
as  the  expected  value  o-f  the  noise  process 
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R  =  E<  vv'  >  (3. 12) 

Then  what  is  needed  are  estimates  o-f  the  measurement  noises 
that  would  be  encountered  in  an  actual  ship  hull.  Using 
studies  o-f  ship  vibration,  the  noise  covariance  .  -for 
displacement  was  taken  as 

r<i,i)  =  3.8  E  -4  meters  squared 

at  amidships.  According  to  Lewis  C27],  the  amplitude 
amidships  will  be  equal  to  about  8.25  the  value  at  the  end 
■for  machinery  induced  vibration.  For  velocity  o-f 
displacement  a  covariance  o-f 

r< i , i )  =  6.8  E  -1  meters  squared  /  second  squared 

at  amidships  was  used.  All  o-f-f  diagonal  values  were  assumed 
to  be  equal  to  zero,  implying  that  the  measurement  points 
are  isolated  -from  each  other. 

The  -filter  was  calculated  using  the  models  -for  the  sea 
spectra,  -forces  and  ship  in  cascade.  This  allowed  an 
investigation  to  determine  the  required  measurements  to 
insure  the  observability  o-f  the  system.  In  addition,  the 
■filter  pole  locations  were  determined.  For  the  Great  Lakes 
Ore  Carrier  Stewart  J.  Cor t ,  described  in  Appendix  A,  the 
poles  were  all  located  at  a  -frequency  o-f  less  than  1.35 
radians  per  second.  The  pole  locations  are  shown  in  Table 
3.1  -for  the  output  matrix 
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c  = 


19  9  9  8  9  8  9 

9  19  8  8  8  8  8 

8  9  19  9  9  9  9 

C  9  9  9  1  9  9  9  9 


(3.13) 


This  output  matrix  provides  measurements  o-f  the 
displacement.  The  resulting  displacement  estimation  errors 
are  given  in  Table  3.2  and  Figures  3.2  and  3.3  show  the 
estimated  value  compared  with  the  simulated  value.  The 
noisy  sional  is  shown  in  -fiqure  3.4. 
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Table  3.1 
Kalman  Fi  1  ten  pol  es 


6.20E-4 

+ 

i     2.07 

7.88E-3 

+ 

1  .08 

9.80E-4 

+ 

7.38E-1 

6.97E-3 

+ 

6.46E-1 

1 .58E-2 

+ 

4.75E-1 

6.39E-1 

+ 

6.50E-1 

5.48E-1 

+ 

6.22E-1 

6.32E-1 

+ 

5.51E-1 

4.75E-1 

+ 

1 .33E-2 

7.37E-1 

+ 

4.92E-4 

6.46E-1 

+ 

4.68E-3 

5.95E-3 

+ 

4.30E-2 

6.08E-4 

+ 

4.SSE-2 

Table  3.2 
RMS  Estimation  errors. 
(meters) 


9.56E-2 
9.36E-2 
1 .54E-1 
9.50E-2 


Prediction    o-f    the    -future    values    o-f    the       states 


qi  yen 


the  current  measured  or  estimated  values  can  be  done  by  the 
equat  i  on 

dxp  =  A  *  xp  (3.14) 

dt 

xp<0)  =  xhat<0) 

I -f  A  is  stable,  the  predicted  value  will  go  to  zero  as  time 
gets  large.  The  objective  then  is  to  determine  how  -far  in 
the  -future  the  prediction  is  acceptable.  I-f  dt  is  the 
prediction  time  step,  the  actual  prediction  becomes 
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xp<t)  =  exp<  A  *dt  >  *  xhat<t>  (3.15) 

which  involves  matrix-vector  multiplication  only. 

For  this  system,  a  number  o-f  prediction  models  were 
tried,  but  none  gave  usable  predictions.  The  -first  model 
tried  was  to  use  the  unstabilized  system,  but  as  expected, 
the  output  immediately  started  growing. 

Next,  the  system  model  was  stabilized  by  using  a  very 
high  measurement  noise  to  calculate  the  -feedback  gain.  For 
this  model,  the  predicted  de-flections  showed  an  initial 
change  in  the  same  direction  as  the  true  value  but  rapidly 
decreased  in  amplitude  to  about  one  tenth  o-f  its  amplitude. 
The  phase  o-f  the  prediction  was  very  close  to  that  o-f  the 
true  de-f  1  ec  t  i  on  . 

Another  method  tried  was  to  stabilize  the  model  o-f  the 
sea  surface  only.  This  resulted  in  a  stable  total  model  but 
no  improvement  was  seen  in  the  prediction.  Here  the 
predicted  value  rapidly  diverged  -from  the  true  value  a-fter 
about  -four  time  steps  and  the  phase  showed  no  apparent 
relation  to  that  o-f  the  true  phase. 

Finally,  a  prediction  was   tried   by   using   only   the 

model   o-f   the   ship.   This  was  even  less  success-ful  as  the 

prediction   immediately   diverged   from   the   true,    o-ften 

increasing    when    the    true   value   was   decreasing   and 

v  i  ce-versa . 

Though  prediction  was  not  acheived  by  any  o-f  these 
models,  this  does  not  mean  that  prediction  can  not  be  done. 
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It  is  one  o-f  the  underlying  ideas  o-f  the  state  space 
technique  that  all  necessary  in-formation  about  a  system  is 
included  in  the  state  o-f  the  system  at  any  time.  Prediction 
models  that  bare  investigation  are:  (1)  stabilize  the  sea 
sur-face  model  by  treating  i  t  as  a  regulator  problem  and 
selectively  choosing  the  elements  in  the  control  weighting 
matrix  to  limit  the  wave  amplitude,  <2)  try  the  same  method 
on  the  total  model,  <3)  a  useable,  i  -f  short  prediction  may 
result  -from  the  second  method  tried  i  -f  the  time  step  is 
reduced,  <4)  use  the  -first  method  tried  but  scale  the 
predi  c  ted  output . 
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CHAPTER  4;  DISCUSS I  ON 

The  problem  addressed  in  this  thesis  is  the  estimation 
o-f  ship  hull  vibrations  -from  a  minimum  number  o-f  measured 
parameters.  These  estimations  are  then  used  to  predict  the 
■future  de-flections  o-f  the  ship  hull. 

The  -first  step  was  to  develop  a  mathematical  model  o-f 
the  ship-sea  system  in  state-space  -form.  This  involved 
three  seperate  models.  The  -first,  which  models  the  sea 
power  spectrum,  was  taken  -from  the  literature.  This  model 
is  based  on  the  ITTC  ocean  spectrum  that  uses  as  parameters 
the  modal  wave  -frequency  and  the  signi-ficant  wave  height. 
In  this  spectra,  the  modal  -frequency  can  be  taken  as 
4.47/Tave,  where  Tave  is  the  average  period  measured  at  the 
calm  water  level.  The  signi-ficant  wave  height  may  be  taken 
•from  table  4.1  i -f  no  measurements  are  available. 


Table  4.1 

wind  speed  (knots)               sig.  wave  height  (m) 

28  371   = 

39  5.1 

48  8.5 

58  11.8 

68  14.6 


The  model  to  include  the  dispersion  relation  -for  sea 
waves  was  developed  -from  a  velocity  potential 
representation  o-f  the  sea,  resulting  in  an  unstaDlc  model. 
The   instability  was   due   to   the   increasing  exponential 
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solution  to  the  resulting  fourth  order  in  time  di  -f  f  eren  t  i  al 
equat  i  on . 

The  third  model  models  the  ship  hull  itsel-f.  This  was 
derived  by  assuming  that  the  ship  could  be  represented  as  a 
simple  beam.  Then  a  composite,  continuous  time,  discrete 
space,  -finite  di-f-ference  representation  o-f  the  simple  beam 
equation  was  used  -for  the  -final  model.  The  model  was 
developed  using  a  second  order  -finite  di-f-ference 
approximation  but  -for  an  actual  application,  a  -fourth  order 
approximation  should  improve  the  accuracy  to  a  significant 
degree . 

Gerritsma  and  Beukelman's  relation  relating  the  wave 
height  to  -force  on  the  ship  was  used  to  couple  the  sea 
model  with  the  hull  model.  Since  in  a  state-space 
-formulation  the  time  derivatives  o-f  the  wave  height  and  the 
de-flection  are  available  in  the  state  vector,  implementing 
the  -force  model  simply  required  additional  terms  in  the 
output  matrices  o-f  the  sea  spectra  and  sea  sur-face  models, 
and  added  terms  in  the  hull  system  matrix. 

Using  these  models,  CTRLC  [263,  a  so-ftware  package  -for 
control  system  design,  was  used  on  a  Digital  YAX-1 1/782  at 
the  Joint  Computer  Facility  to  do  the  simulations,  develop 
the  -filter  gains  and  evaluate  the  sensitivity  to  errors  in 
the  model.  The  estimation  o-f  the  de-flections  o-f  the  hull 
could  be  done  well  with  the  state  space  model  developed. 
The  prediction  o-f  the  de-flection  to  some  time  in  the  -future 
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was   not   successful   using   the  models  developed  here,  but 
with  improved  models  it  could  be  possible. 
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APPENDIX 
PREDICTION  RESULTS 

The  several  methods  o-f  prediction  attempted  gave  less 
than  satisfactory  results  when  compared  to  the  simulated 
signal.  Figure  A.l  shows  the  prediction  that  resulted  -from 
using  the  system  model  stabilized  by  using  a  very  large 
measurement  noise.  This  resulted  in  a  predicted  value  that 
was  much  larger  than  the  simulated  value.  The  same  general 
pattern  resulted  when  the  measurement  noise  was  varied  -from 
a  low  equal  to  the  expected  measurement  noise  to  as  high  as 
could  be  computed.  Nor  did  changing  the  starting  point  o-f 
the  prediction  with  relation  to  the  phase  of  the  simulated 
signal  improve  the  prediction. 

The  next  -figures  show  the  prediction  resulting  -from  a 
model  constructed  by  stabilizing  the  sea-sur-face  model 
be-fore  -forming  the  total  model  .  In  Figure  A. 2,  the 
predicted  signal  is  very  similar  to  the  prediction  in 
Figure  A.l.  I-f  just  the  -first  -few  time  steps  are  plotted, 
it  can  be  seen  that  the  predicted  signal  does  not  track  the 
simulated  signal  even  -for  a  short  period  o-f  time,  Figure 
A. 3. 

Finally,  Figure  A. 4  shows  the  results  o-f  using  only 
the  ship  model  to  do  the  prediction.  As  can  be  seen,  this 
prediction  has  no  apparent  relation  to  the  simulated 
si  gnal  . 
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F  i  gure    A. 3 
Stabilized    sea-surtace,    -first    9    time    steps    o-f    prediction 
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F  i  gure  A. 4 
Prediction  using  only  the  ship  model 
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